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INTRODUCTION 


This brief report will review the progress made under the subject Grant in the 
period 11/1/92 - 5/31/93. The research involves the continued development of the 
Large Angle Magnetic Suspension Test Fixture (LAMSTF) and also the 
recommissioning of an additional piece of exisiting hardware. 


REVIEW OF CURRENT AND ONGOING WORK 

During the period in question, the initial configuration of LAMSTF was 
completed and made routinely and reliably operational. A digital phase advance 
controller was completed and documented. The goal of a controlled 360° rotation was 
been achieved. Work started on the recommissioning of the Annular Suspension and 
Pointing System (ASPS). A more detailed breakdown of work completed follows : 

(i) Modelling. A series of simplified eddy current measurements have been made 
which show promise in resolving the related dynamic modelling problem. A new 
version of the computer code GFUN has been supplied and has been installed and 
tested on the SCB SUN network. 

(ii) Position sensing. Design work was started on a revised optical system, 
intended to overcome continuing difficulties is establishing and maintaining calibration 
of the existing system and also to incorporate roll sensing. 

(iii) Controller. The digital controller previously developed under this Grant is 
presently based on 486-class PC hardware, implementing classical dual phase-advance 
compensators. Preliminary documentation has been provided. 

(iv) Support of the Second International Symposium on Magnetic Suspension 
Technology. Calls for papers and exhibits were prepared and distributed. Submitted 
Abstracts were reviewed and author kits sent. 

(v) Recommissioning of the Annular Suspension and Pointing System. The 
hardware has been uncrated, inventoried, positioned in a laboratory and reviewed for 
condition and ease of recommissioning. Studies of the modifications required to permit 
suspension in a full gravity field and a review of the control problems have been 
completed. Copies of the student project papers generated are attached as Appendices 
to this report. It was concluded that simple shims added to the pole faces would greatly 
increase the force capability and easily permit suspension in a gravity field. The 


increase the force capability and easily permit suspension in a gravity field. The 
fabrication of the pole shims required to permit such suspension is underway. 


PROBLEM AREAS 

GFUN has continued to be a rather troublesome code, although the SUN version 
does execute considerably faster than the previous VAX version and has much superior 
graphics. The current graduate research assistant, Lucas Foster, is becoming proficient 
with the code. 


PUBLICATIONS DURING THE CURRENT GRANT PERIOD 

Groom, N.J.; Britcher, C.P.: A Description of a Laboratory Model Magnetic 
Suspension System with a Large Angular Capability. 1st IEEE Conference on Control 
Applications. Dayton, OH, September 1992. 

Groom, N.J.; Britcher, C.P.: Stability Considerations for Magnetic Suspension 
Systems with Electromagnets Mounted in a Planar Array. NASA TP-3229, December 
1992. 

Ghofrani, M.: Approaches to Control of the Large Angle Magnetic Suspension 

Test Fixture. NASA CR-191890, December 1992. 
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TNTRODUCTION 


Levitation, as far as this report is concerned, is the act of 
holding- up an object with no visible support by means of 
electromagnetic suspension techniques. Since the 1930 's there has 
been research in single axis suspensions utilizing a magnetic 
bearing station as shown in Figure 1. 



SINGLE MAGNET SUSPENSION 


FIGURE 1. 

This research has blossomed into the idea of multi-axis 
suspensions. Multi-axis suspension has several advantages over 
single axis system, in that it provides control of an object with 
precision in two or three orthogonal axes. In this report, we 
discuss the primary use of magnetic -bearing suspension and it's 
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rel evance to what was formally known as NASA's Annular Suspension 
an d Pointing System (ASPS) . This system is an experimental 
pointing system with applications for the space shuttle and the 
space station programs. 

The objectives behind this magnetic suspension research 
project are to provide insight to the use of the ASPS 
configuration, to control the solar panels of the space station. 
This is important to maintain the correct position of the panels in 
relation to the sun and orbiting space station for the continuous 
supply of solar energy. Since the panels are suspended, they can 
be aligned with minimum outside interference . The approach of 
using magnetic suspension technology guarantees mechanical 
isolation since there are no contacting surfaces. This isolation 
reduces vibration transmission and mechanical wear which in turn 
extends the life of the payload and of the carrier. It should be 
noted that ASPS has a high pointing accuracy along the line of 0.01 
arc -second. 

This research will be done in a laboratory setting by 
incorporating five bearing stations and one motion control station 
(See Figure 2) . We will attempt to suspend an object of dead weight 
similar to that of a solar panel. The long term applications may 
include deep- space navigation, fire control in weapon systems, and 
an improved mass transit system. 
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THEORY 


The principle behind electromagnetic suspension in this 
project is simply explained by magnetic flux theory. The magnetic 
flux through any closed figure is the product of the area by the 
average component of magnetic induction normal to the area. Flux 
($) is defined by : 

<P = k f B dS 

J 3 


where k is the constant of proportionality between field and flux 
density (B) . This application can be expanded to represent the 
principles surrounding single axis suspension. An object is 
supported, against the force of gravity, by an electromagnet in 
which the current is controlled electronically in response to a 
position signal. From Earnshaw' s Theorem, it is known that the 
stability of the system is dependent upon the feedback of this 
position signal because without control of the current it is 
fundamentally unstable . To achieve stable suspension, it is 
necessary to regulate the current in the electromagnet using 
position feedback of the object to be suspended. 

The force of attraction between two objects is given by the 
formula, F m = B 2 + 2\i q * Area, where B represents the magnetic flux 
density, [i 0 is permeability of free space, and the area is 
represented by the cross-sectional area of one pole. This formula 
comes from the derivation of Maxwell's stress formula. The 
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electromagnets used in suspension systems are under the influence 
of an -air gap". The air gap denotes a gap left in the magnetic 
material; a abort gap is usually left in the core material to 
prevent saturation of the magnetic material by the dc current. 
This implies the gap flux density is directly proportional to the 
ampere turns -MI" and inversely proportional to the length "g" , gap 
length. The formula then becomes, F m - (l/2p 0 ) * <Uo * MI/g; 

Area, as seen by the following derivation. 

’ ZL 


where H m is the magnetic field 
material . 

Maxwell's equation states: 
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The governing equation states 
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However, in this system we have A = 2* (width* length ) , as with 
the air gap, taking into account two pole faces. 


F „ = Vo 


(NI) 2 

9 2 


Area 


The implication of this gap distance is dependent upon 
several separate factors, the weight of the suspended object, the 
induced current, and number of turns on the coil. (See Figure 3) 
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Figure 3 

With a 

single axis suspension system it is obvious to see that there 
only one position on the vertical axis at which the magnetic lift 
force is equal to the weight of the object. Any deviation to this 

position will result in the object's displacement. 

The translation motion of a free body in space will have three 
degrees of freedom associated with it as well as three degrees 
freedom relating to rotational motion. Again, based upon the 
classical force equations, the motion equations can be derived for 
this multi -bearing system. The motion of the system can then be 
described by six non-linear second order differential equations. 
The non-lineararity could be due to bearing, characteristics and/or 
rotations due to angular momentum. The equations can be linearised 



after utilizing Maxwell's stress formula to find the magnetic 
stress tensor. The resulting matrix of the tensor provides values 
for all components of magnetic stress along each of the coordinate 
axes. The rigid body force acting on an object is obtained by 
integrating these components over its bounding surface. As far as 
the rotation due to angular momentum is concerned, a summation of 
the moments is calculated and then linearized. It should be noted 
that for efficient operation of suspension systems utilizing dc 
magnets, small air gaps which can be precisely controlled are 
required. 

In complicated systems such as the ASPS, multiple bearing 
stations are utilized. There are several difficulties associated 
with multi -magnetic systems to include the management of six 
separate air-gaps . The integrated control portion will be 
addressed in a separate report. 
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DESIGN APPROACH 


Using the related force equation, F = /i 0 + [(NI) 2 + *Area, 
we are to modify the existing multi-bearing system to meet the new 
requirements assigned. From the design documentation, the weight 
of the rotor was given as 47.6 lbs or approximately 211. 7N. This 
implies that each bearing station should support 16 lbs or 7 1.2N. 
However, to allow for errors and possible addition of payload 
weight, a 1.3 safety factor is incorporated. Each station should 
hold 20.63 lbs or a force of 91.75N. An additional factor of 1.5 
is included, for control purposes , to give a sum total of 3 0.94 lbs 
or 137.63N per axial station. This 1.5 factor is the peak force 
prior to saturation of the magnetic material. Since the axial 
magnetic bearing station can hold twice the amount as the radial 
bearing station, we chose to redesign the axial station first in 
order to get a reasonable approach to the new requirements . In the 
formula above, there are five variables that can be changed. 
Systematically we selected the variable to be altered and made the 
appropriate calculations . 

In the first design, we regarded the (NI) 2 term to be a single 
factor and the only variable to be considered while all the others 
were held constant from the original design specifications . This 
calculation indicated that the "NI" term had to be increased by 
54.5 % to meet the required specifications. The "NI" term can be 
brok en into it's two components, one involving current and the 
other number of turns. The number of turns would be very difficult 


8 



to increase due to the replacement of the entire coil and the 
protective coating on the wire. Since it would not be cost 
effective to increase the number of turns, the alternative is to 
increase the current. This also is not feasible because increasing 
the current will cause other problems in the system. These include 
overheating of the bearings, melting of the wire, and magnetic 
material saturation . The comparison of the original and new design 
specifications are shown in Table 1 and Graph 1. The only 
advantage of this approach would be that the force would increase 
substantially since F <x (NI) 2 . 

In the second design, we regarded the area to be a single 
factor. To meet the new specifications , the area was calculated to 
be 1.4E-3 square meters and an increase of 79.3 % as shown in Table 
1 and Graph 2. We took into consideration an increase in length, 
width, and a combination of both dimensions. This approach was 
again not feasible due to cost constraints. It would be very 
expensive to rebuild the bearing station magnet to meet the new 
specif i cations . 

The third design took into consideration the gap width. The 
original specifications set the nominal gap width to be 
approximately 6.35mm and at the time of this research this was an 
acceptable distance. Since that time, small-air gap theory, 
commonly used for magnetic bearing for machinery, has decreased the 
gap distance to 1/100 of an inch. This new concept led us to the 
conclusion that the movement of the bearing to decrease the gap 
width would be the most cost effective and easiest design change. 
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Taking in to consideration the original design specs and keeping 
all variables constant , we calculated the gap width to be 
approximately 3.466mm which is a 83 * decrease in distance (see 
Table 1 and Graph 3) . If this decrease should cause saturation or 
over heating problems, the other design variables will have to be 
altered. In the testing portion of this report, it is shown that 
there was no overheating or saturation due to the change in the gap 
width . 

The radial bearing station designs will be minimally modified 
because the translational motion is negligible . From the original 
design, we determined that the variables could remain constant with 
the same change in gap width as the axial bearings. 


TABLE 1 - 



ORIGINAL 

DESIGN 

* CHANGE 

N--# OF TURNS 

2148 

4720 

+54.5 

I- -CURRENT 

1.4 

2.2 

+36.4 

A-- AREA 

2 . 9 0X1 0 " 4 vc? 

1.40X10 ' 3 vc? 

+ 79.3 

g- -GAP WIDTH 

6.35mm 

3 . 47mm 

-83.0 
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TESTING AND DESIGN EVALUATION 


The original ASPS design and testing was accomplished 
in the mid 1970' s to the early 1980' s. Since that time, the ASPS 
program has been on hold. The purpose of this experimentation is 
to verify the operational capability of the rotor assembly under 
new design constraints . 

In the new design, the gap width was decreased to increase 
the lift capabilities . As mentioned in the design section, 
saturation of the material and overheating of the coils could be 
two problems encountered by decreasing the gap width. A 
mathematical exercise can be used to prove the material does not 
saturate . Knowing that the flux is a constant through out the 
circuit, flux density in the smallest cross-sectional area can be 
calculated . The following calculations will show that saturation 
should not be a problem with the new design. 

From previous derivation: 

B a = ~~r ~ = 0.389 Tesla 
a 2g 

<P = constant = B a A a = B , A , 

« d J. ± 



A± will be the smallest cross-sectional area, which will yield 
the largest "B"; this will show that saturation does not occur. 
A a = 1.55E-3 = 7.76E-4 B a = .389T 
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B i = - Hfg I (0.389T) = 0.78T 
7 . 76E-4 

As shown in Figure 4, B sat for the material is given as 1.6T. 
The calculated as ,4T and the "B n for the smallest 

cross-sectional area is .8T, which is well below the B level . 
There should be no saturation problems due to the new gap width. 



56 112 224 

AMPS/METER 


The second consideration mentioned was the overheating of 
the coils. This was proven in the lab by running a current of 
1.5 amps through the coils for twenty minutes. A thermocoupler 
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was connected to the outside casing of the coil and a piece of 
foam rubber was positioned around both to insulate them from the 
outside air. As seen in Graph 4, the increase in temperature was 
minimal. The current was increased to 2 amps and smother twenty 
minute test of time versus temperature was performed. Graph 5 
indicates the increase in temperature was again minimal. 

Each of the bearing stations were tested for their magnetic 
force by using a 1.5 amp current and a magnetic sensing devise. 
This proves the magnets are in operating condition . 

From these three simple experiments, we believe that our 
design will function as desired under the specific design 
criteria. 

The actual testing of the rotor assembly was accomplished by 
the following steps: 

1) The schematics indicated the power to the axial bearing 
stations came from plug J402. Each bearing station is 
designated by either A, B, or C with pins 1, 9, and 20 being 
the input and 2, 10, and 21 being the output. To have each 
bearing station operating simultaneously, the output of A- 
pin 2 was connected to the input of B-pin 9 and B-pin 10 was 
connected to C-pin 20. A-pin 1 was connected to the power 
supply and C-pin 21 was connected to the return. 

2) To obtain the appropriate current, two power supplies 
were connected in series. It should be noted that with them 
connected in series, the current was monitored by an 
external amp meter to insure the exact current. We adjusted 
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the current to approximately 1.6 amps and manually lifted 
the rotor until it was magnetically held by the top portion 
of the bearing assembly. We then decreased the current at 
0.01 increments until the rotor released itself at 
approximately 1.35 amps. 

3) The desired air gap was maintained by inserting a 
non-magnetic material (aluminum) shim to simulate an 
air gap of 3.46 millimeter. Therefore, when the rotor 
was raised, the assembly was operating at the designed 
criteria. The results are shown in the accompanying 
ph o t ograph s . 
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CONCLUSION 


Based on the design calculations, the assembly operated 
within acceptable limits. The current designed for was 1.4 amps, 
however, the assembly operated correctly at 1.35 amps. This 
indicates that the assembly could operate at lower currents and 
also implies there would be no problem with overheating of the 
coils or saturation of the magnetic material: 

The follow-up research for the magnetic bearing system 
should include integration of the roll -motor assembly with the 
findings of this report. This would give control of the entire 
bearing system to roll -motor assembly instead of the few pins 
referenced in the experimental section of this report. A 
computer program could be used in conjunction with the control 
motor for ultimate control of the system. 
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List of Symbols 


A: cross-sectional area of one side of the actuator 

a: assigned constant 

B: magnetic flux density in the gap 

Bm: magnetic flux density inside the actuator 

b: assigned constant 

c 2 : assigned constant 

F(x, i) : attractive force between the actuator and the rim(suspended 
mass) 

G p (s) : plant dynamic equation or open loop transfer function 

g, , g 2 , g: gap distance 

g„: equilibrium distance 

I : input current of the actuator 

lb/ I.: equivalent bias current 

i: controlled current 

L: inductance of the actuator 

L c : inductance of the coil 

Lj: inductance of the gap 

1: the length of actuator 

m: mass of the rim (suspended mass) 

N: number of turns of he coil 
R: actuator resistance 
s: Laplace transform variable 
V: input voltage of the actuator 
V b : equivalent bias voltage 
v: controlled voltage 

x: controlled distance or small perturbation distance from 



equilibrium 


H : permeability of the actuator material 

K : permeability of free space 

<t> : magnetic flux 
L.H.P : left half plane 

ASPS : Annular Suspension and Pointing System 



I. Introduction 


A. Problem Definition 

Frictionless electromagnetic suspension and levitation has 
attracted much attention since 1970. Applications include high- 
speed machine tool spindles, ultra-centrifuges, high vacuum pumps, 
and fly-wheels for energy storage. Methods of producing 
electromagnetic suspension and levitation include controlled DC 
electromagnets, diamagnetic materials, superconductors, hybrid 
systems, and tuned LCR circuits. A comprehensive review lecture of 
electromagnetic suspension and levitation techniques can be found 
in reference 1. 

The technique of suspension and levitation with controlled DC 
electromagnets is the most advanced and successful at this time. 
Many investigations are underway worldwide. Advanced ground 
transportation schemes, contactless bearings for ultra— high speed, 
and gyroscopes have been successfully demonstrated by many groups 
of researchers. 1 

The Annular Suspension and Pointing System (ASPS) developed by 
the Flight System division of Sperry Corporation 1 is a six-degree 
of freedom payload pointing system designed for use with the space 
shuttle. This magnetic suspension and pointing system provides 
precise controlled pointing in six— degrees of freedom, isolation of 
payload-carrier disturbances, and end mount controlled pointing. 
Those are great advantages over the traditional mechanical joints 
for space applications. More detail discussions of the magnetic 
suspension joints and mechanical joints can be found in reference 
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2 . 


Figure 1 and 2 show the ASPS designed by Sperry Corporation. 
It consists of six actuators, three for vertical movements, two for 
radial movements, and one for tangential movements. By the coupling 
and decoupling matrices (figure 3) 2 , we can carefully decompose the 
command signal of each degree of freedom to each actuator 
individually. In other words, the coupling and decoupling matrices 
change the six-degree of freedom ASPS control system to six single- 
degree of freedom ALPS control systems. Hence, we can design each 
control loop separately. 


figure 1. 





figure 3^. 
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B. Project Objectives 


(i) Understanding the existing ASPS in the lab. 

(ii) Model the dynamics of a single DC controlled ALPS 
actuator as accurately as possible. 

(iii) Re-design a controller for the single degree of freedom 
ALPS control system to achieve the highest stiffness as possible. 
{Highest stiffness will have the lowest motion in response to 
external forces} 

C . Summary 

In this design, we first analyzed the assumed model of the 
single degree ASPS bearing actuator, and obtained the plant 
dynamics equations. By linearizing the plant dynamics equations, we 
designed the cascade and feedback compensators such that a stable 
and satisfied result was obtained. The specified feedback 
compensator was computer simulated with the nonlinearized plant 
dynamics equations. The results indicated that an unstable output 
occurred. In other words, the designed feedback compensator is 
fail. The failure of the design is due to the Taylor's series 
expansion does not converge. 
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n. Modelling of the single degree ASPS bearing actuator 

A. Assumptions of the model 

The single degree ASPS bearing actuator consists of two pairs 
of magnetic coil elements, mounted in opposition, to control the 
rim (suspended mass) along a single axis. Figure 4 shows the 
configuration of the actuator. The magnetic coil elements have 
current biasing superimposed by a controlled DC voltage source to 
produce a force to suspend and point the rim. For fine pointing 
application, the gap distances between the rim and actuators are 
kept to a predetermined value (<0.3") . 

figure 4. 
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We derived the equation for the single ALPS bearing actuator 

with the following assumptions: 

(i) The force of attraction between! magnetized bodies is given 

by F = (l/2u) * B 2 * area. 

(ii) The magnetic flux density is uniform between the gap, or 

gap distance « size of actuator. 

(iii) The rim is a perfect conductor. That is, it does not 
support any magnetic field strength (H) inside the rim. 

(iv) There is no coupling effect between each coil elements. 

(v) The controlled electromagnet behaves linearly , and there 
is no loss or fringing of magnetic field. 

(vi) The rim has a plane area over the magnetic coil 

assemblies. 

r. Derivation of ASPS Dynamics 
a) B-field in the air-gap. 

From magnetic circuit theory, the total current linked by the 
path of an N-turn coil is given by 


Assume 

(i) the magnetic material is approximated by B m = Mm^m + 
B c where B c is a constant. 

(ii) gap distant at path (1) = gap distance at path (2) 
ie. g x = g 2 =g 

(iii) from boundary condition, B inside the coil ~ B in 
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the gap. 


Therefore, equation (l) becomes. 


+ H 2&2 + H m l = Ni 

B - B (because the 

** 2 Hg + —2 °i = Ni rim is a perfect 

conductor) 

- 2Hg + £°1h = Ni + 


As the term B 0 (l//x m ) can be equivalently assigned to 
a magnetomotive force (mmf) (B 0 //x m ) 1 = NI 


So 


2 Hg + = jvi + wj o 


Thus the magnetic flux density at the gap is given by 


B = \i a H = + Jp) 

2g + 

As for the particular material we used for actuator, n m »n 0 . 

Hence 

B • MU 1 J °> 
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b) Relationship between F, I and g 

By considering the stored magnetic energy, Bohr 3 and Hayt 6 were 
able to relate the magnetic attraction force to the magnetic flux 
density and cross-sectional, 

F. = . area 

The geometry we used is similar to Humphris 7 and Groom*, figure 
4, which have two electromagnets positioned opposite the rim. This 
kind of configuration is more linear if we separate the magnetic 
flux density into the controlled and bias components. 7 

Assume no coupling effect between the two actuators, by 
equation (2) : 

1 2 ^ ' = 2 ^ 

Let = N 2 = N, Ij = l 2 = i o , ij = -i 2 = i = controlled current. 
Therefore, the total force acting on the rim is given by 
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= Ft - F 2 - J24L (B i - Bi > - < 8 , * B,) (B, - B a ) 


2|» 


J„ ♦ i I„ - i. 

Consider B 1 - B 2 = — 


(gr a ♦ x) (J c + i) - (go - *) (J 0 - i) 
2 grj - x 2 


] 


2 (go - x 2 ) 


^ — [g 0 (J 0 + i - J„ + i) + *(*<, + i+ I 0 ~ i >3 

«*i ♦ «.* - M<&* - 

2 (gl - x 2 ) 9o ~ x 2 


Similiarly, B 1 + B 2 


* s . , MtlLLi i ±S±i - —^— („A * *i> 

2 g 0 ~ x g 0 + x g\ - x 2 


a . . A r Vjf(9j + I a X ) -, r V‘o N ^o I 0 + **> 

Thus F = — (B x + B a ) (B x - B 2 ) = — [ j r ] t j “ 

H 0 Vo gl - x 2 gl - x 2 


iij^MgJ + J^) (gr 0 J 0 + *i) 
(<7 2 -* 2 ) 2 


By the Taylor's Series Expansion at the equilibrium point (x 0/ i 0 ) , 
we get 
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F(x, i) 

dF 

as u 


F[x ol i 0 ) + <x - x 0 , i - i 0 )*VF| Uo> ig) 

7 + xi > + + *«*) 1 

(^O -X 2 ) 2 


= T-T 2 — 777 I Iffi - x 2 ) 2 [I 0 {g 0 I 0 + xi) + i(gr c i + J^x) ] 
ox (9o ~ x 2 )* 

~ 2 (g 2 - x 2 ) (-2x) ( g 0 I 0 + xi) (gr D i + 1^) ] 
at equilibrium point ( x 0 , i 0 ) = (0, 0), 

F(.x ot i 0 ) = 0 

dF, = 

dx tJf °' io> 


3F, = 

di (Jf °' ie> ^2 


So 


F( X , i) - X. * i * ij 






2 


(3) 


c) V— I relationship of the actuator. 
Recall equation (2) : 


B _ |i jHi +_ J 0 ) 

2g 
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By definition, *=NBA and 


L=(d*/di) 5 


Therefore, the 


inductance of the gap is 





N^Ay 0 

2g 


and the inductance of the whole circuit is 


L = L t + Lj + L c = 


N 2 A\i a j N 2 A\i 0 ^ r 

■■ + ” T j-t — 

2 g x 2g 2 

n g^ogo + L 


(4) 


(by 4) 


By Kirchhoff's voltage law, we have 


V = Ri + (Li) 

at 


Ri* 

dt dx dt 


We previously separated the voltage, current and gap distance 
into the bias components and controlled components. 

That is Let V = V b + v 

i = I b + i 

and x = g 0 + x 

Therefore, 
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V„ ♦ V = R(l b ♦ i) * I*i * * l". *E 

* * dfc dx dt b dx dt 


As, at equilibrium position, V b = RI b , assume i«0, x«0 
Therefore, 


• dL dx « 
2 -==*-== * 0 
dx dt 


so v = Ri + 

dt b dx dt 


_ . . rfr 2.W 2 .ALl-(7- 

Consider -== = 

S^o - x 2 


dL^dx T 2^ 2 A|i 0 5r 0 


dx 


T . V 4 J . i UnAi I riw ri ^ 

^ S St - z “—- * X *- T - ■ 0 


(?!-x 2 )* dt 


Thus v(t) - Ri ( t) + L-^-i(t) 

dt 


taking the Laplace transform on both sides, we get 


V(s) = RI (s) + Lsl(s) 

or J(s) = — 1 V(s) 

R + Ls 


Recall equation (3) 
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F(x, i) 


= A»j£ll x + A^I 0 ± 


9 l 


9o 


By Newton Second Law F(x, i) = m 


d 2 x 

dt 2 


so. 


dt‘ 


3 

9o 


9l 


Taking the Laplace Transform on both sides, we get 


ns’Xis) - * A ^‘ T ° 1 


9l 


9o R + Ls 


V(s) 


or 


X(s) _ 

via ) 


9o 


(R + Ls) (ms 2 - 

9 l 


( 6 ) 


AV pH 2 Ip 
9 pfnL 


( a + -) ( s 2 - ) 

L n. < 7 0 3 


Let a - 


9 lmL 


C 2 . 


2 T 2 
o 


m g\ 


Therefore the plant dynamics of the ALPS actuator are 


,0 4 m - 


Vis) ( s + b) ( s 2 - c 2 ) 


Which is similar to the plant dynamic equation obtained by 
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Kilgore 9 and Jay a want 3 . 

Referring to Groom 2 , the values of those parameters 


Io 

0.57 Amps 

A 

1. 1400918*10‘ 3 m 2 

V-o 

4 n *10' 7 H m' 1 

N 

1386 turns per coil 

m 

7.19712 kg 

g« 

0.00762 m 

R 

8.0 Q 

L 

0.1805899 h at g„ 


.-. a * 20.79 
b = 44 .3 
C 2 = 280.8 

Thus, the open loop transfer function is 

20.79 

Gp s 3 + 44.3s 2 - 280.8S - 12439.59 


are. 
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d) Discussion of the plant dynamic of ASPS actuator 
The open loop transfer function is a third order, type zero, 
all poles plant system. The characteristic equation also contains 
one positive real root, so this plant is not BIBO stable. The pole 
zero diagram of the plant is shown in fig 5. In order to move the 
open loop unstable root into the stable region, we need to add 
zeros in the left half plane so that the locus are pulled into the 
stable region. In other word, a reshaping of the root locus 
(compensator) is necessary. 



figure 5 
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IH. Compensators Design 


A. Cascade Com pensator 

The approach we used in the cascade compensator is to achieve 
the goal ( shift the root locus to the L.H.P. ) with minimum 
complexity. For the simplest case, by the knowledge that a pole 
will pull the root locus to the right, a zero will pull the root 
locus to the left and a pole-zero pair close to the origin will 
decrease the steady state step error (G p is type zero) , we tried the 
general lag— lead cascade compensator , 


G C (S) 


(s + 


4 ) 


= A 


(s 


a T, 




with a =10, and gain A. 10 The lag component was fixed at 

(s+ 0.05) / (s+ 0.005) and the lead component was moved along the 
real axis. Some root locus results are shown in appendix I. After 
studying the results, we decided that we needed to increase the 
compensator complexity in order to meet the design specifications. 

Since the lag component only affected the steady state error, 
for simplicity, we tried the dual phase advanced compensator with 

a=10. 


G c (s) 


= ( 


3 + 


S + 



2 


Some root locus results were shown in the appendix II. 
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A block diagram is shown in figure 7. We selected a damping ratio 
5 = 0.7, and the maximum natural frequency. The figures of merit 






poles: -86.15 ± 87 . 89 j , -93.24, -620.46 

additional gain K = 1.011406 E 6 

steady state step error =0.06 

rise time = 0.009 sec. 

peak overshoot = 1.42 

peak time = 0.025 sec. 

settling time = 0.05 sec. ( for 5% ) 

gain margin = 16.75 dB 

phase margin = 180° 

stiffness, <*> n = 15 kN kg'm 

stability region : 6.0939 E 4 < K < 6.8041 E 6 

The step response was shown in figure 8, and the Bode plot was 

shown in figure 9. Those results are obtained by matlab. 

figure 7 . 
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step response of cascade compensator 



figure 8 
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B. Feedback Compensa tor Design 

The cascade compensator will be converted to a digital 
controller later, and a computer will be involved to control the 
plant. Therefore, it is natural to design the compensator using the 
state feedback technique. This technique is flexible and convenient 
to implement. A brief derivation of the design procedure is shown 
in Appendix III. In the case where some state variables are not 
accessible, an observer (estimator) may be used. Observer design 
procedures are also shown in the Appendix III. This material are 
come D'Azzo 10 . In this design, we used the full state feedback 
technique . 

A state space representation of the open loop plant is shown 
in the figure 10. 


figure 10 

State space representation of the open loop plant 


X 


~y. d 


j. * 
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*1 

Xr 
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o 
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I Xz.i 
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The controllability matrix is, 
Mc=[bjAb|A 2 b] 


0 

1 

» 0 

1 

1 

t 

! 20.79 

1 

1 

0 

1 

! 20.79 
1 
1 

I 

! -920.997 
1 
1 

20.79 

1 

! -920.997 
1 
1 
1 

1 

! 46637.99 

i 


as det( Me) . 0, so Me has full rank. 

The observability matrix is 
Mq= [ C T | A t C t | ( A t ) 2 C T ] 

I I 

i ! o ! o 

I I 

I I 

I 1 

o ! l ! o 

l I 

l l 

i i 

o ! o ! l 

I t 

I t 

I I 


So Mq has full rank. Hence, this system is completely 
controllable and completely observable. 

Motivated by the performance of the cascade compensator, we 
selected the poles of the control ratio to be -86.15 + 87.89j, 

-100, which give us a good step response. A block diagram is shown 
in the figure 11. The figures of merit are, 
gain = 72.8546 E 3 
k t = 1.01 
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0.021561 


k 2 = 
kj = 0.0001505 

poles : -86.15 + 87.89j, -100 
steady state step error = 0 
rise time = 0.03 sec. 

settling time = 0.045 sec. (for 5 %) 
gain margin = 13 dB 
phase margin = 160° 
stiffness = 15 kNkg' 1 m' 1 

A step response and Bode plot were shown in figure 12 and 
figure 13 respectivity . Those results are obtained by matlab. 

figure 11. 
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Amplitude 


step response 



figure 12 
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figure 13 
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IV Computer Simulation 


We simulated the state feedback design by a four order Runge 
Kntta Method in the time domain. The Runge Kntta Method is a 
special version of the general Taylor's series expansion. 

For a general different equation. 


dX 

dt 


= 2( x ) 


Let, 


*!«*<*> 

- 1 ( X + 1 h ) 

- 2 ( X + 1 h t 2 ) 
2 \ - i ( x + h f 3 ) 


The next x can be approximated by 


Xnext = X + | [ F l + 2 ( + *3 ) + ] 


with an error of fifth power term of the Taylor series 
expansion. The h is the increment step side. 

Recall, 


32 



mx= 


A[ i 0 N 3 


(9o I o 


+ x i) (g Q i + 


I a x) 


(g „ 2 - x 3 ) 2 


V 


- Ri + L ik 


Let 


x 1 = x = output 

dx 1 
* 2 = ~dt 



v = u = input 


So 


*1 = *2 


A H 0 N 2 


x 2 = -S’ = 


m 


(g a i 0 + *1 x 3 ) (g a x 3 + J 0 x x ) 


( 9 T 0 2 - X 1 2) 


2 \ 2 


*3 = 


v + — LZ 

T 3 L 


x 4 “ l 


Let the control law be 


u = r - ( Jfx + Jc 2 x 2 + ic 3 x 2 ) 


where r is the reference input. 
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second 


In the computer code, we need to estimate the 
derivative of x, which is achieved by 

( x 2pt*wat “ X 2 previous) / time interval 

For the linearized plant dynamics equations, slightly 

modification of the computer codes can do the job. The results of 

computer simulation are shown in figure 14 to 19. The computer 

program was written in Pascal language, and is shown in Appendix 
IV. 
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V. Discussion of results 

The computer simulation showed that the step response 
increased linearly with time, figure 17. That is, an unstable 
response occurred. 

Consider the plant dynamics equations, recall 


A \i 0 N 2 ( g„ Ip * x i ) ( 9o i + J o x )_ 

X ' 


( 9o -X 2 ) 2 
As for | x | < | 9 0 I 
1 1 


(<T 0 2 - * 2 ) 2 


& O 


( 1 - -^) 2 
9 0 


= ^o ‘ 4 t 1 + + (_ f ^ )2 + " ‘ 32 

y 0 y o 

= „-4 r 1 + 2(-^i) + 3(-^) 2 + higher terms ] 

0 g 0 2 z 2 

(g a l 0 * x i) + i q x ) 

50 ■ * 7 2 \ 2 

(ffo - * > 

= ^o" 4 <^o I o + * i) i + J o X) 

[ l + 2 ( — ) + 3 ( — \ ) 2 + higher terms ] 

9 2 9o 

= g,,-* (g 0 J 0 2 x + g 2 I 0 i + 9 0 * i 2 + T o x 2 2) 

[ 1 + 2(^i) + 3(^i) 2 + 
g- Q sr 0 

= 9 0 -* t ( 9 0 *o 2 x + <7 e 2 

+ (2 i°! x 3 + 3 I Q x 2 i + g 0 x i 2 ) + iiigiier terms ] 
3o 


Obviously, the coefficient magnitude of the third terms in the 
Taylor ' s series expansion is larger than the coefficient magnitude 
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of the first term. Thus, at least, we need to include the third 
terms in the compensator design. However, we cannot use the 
conventional linear design theory in this situation. 

Even thought we include the third terms in the design, the 
coefficient magnitude of the x“ terms in the series expansion 
increase without bounded ( because g c < 1 ) . So, the Taylor's series 
expansion does not exist. 
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VI. Conclusion and recommendation 

The plant dynamics equations are nonlinear, and the 
conventional linearization does not work. We recommend to design 
the compensator without linearization in the time domain. 

As, in general, if we close the loop, 

=?{ x , u ) 
at 

u = g {z , x ) 


define PI to be 

f“ t (x x - r ) 2 dt 

J 0 


and, minimizing the PI which is subjected to the constraint 
equation, 

x , g ( r , x ) ) 

at 

, by the Langrange multiple method. The mathematics is too 
difficult to carry out, and the analysis is left for interested 
reader only. 
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Appendix I 


Root locus result of general lead lag cascade compensator 



— 300 —200 -lOO O *1 DO 200 300 


R*al AxTv 


Co — cC*-*-oas)/(»-»-o.oos)3CC*-^20>/<*H-2so)i 
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Appendix ^11 

State feedback a& siqn procedure 10 


Let the open-loop plant transfer function be 


g p = — 

p s 3 + a 2 s 


v = - 

2 + a[ s + a 0 u 


Let, for standard notation. 


*1 = y> x 2 = -JZ 


_ dy x = d 2 y 


dt ' 


d t‘ 


dx 1 _ dy _ „ dx 2 _ d 2 y = ^2 m 

— rr - —rr ~ x 2 ' ~TZ ~ m x 2' w*- 


dt dt 


dt dt‘ 


dt 


d 3 y 

dt 3 


and 


dX 

dt 


= A X + h u 


y = [ 1 0 0 ] X = £ T X 


Let the control law be 


u = r - k T X 

whexe k T = [ k 2 ic 3 ] 


Consider 
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H(S) 


k^_Xis) m k*X 
X(s) c T X 


k 1 x 1 + k 2 x 2 + Jc 3 x 3 
— — 


H(s) 


k 1 x x + k 2 sx t + Jc 3 s 2 x 1 
_ 

k 2 + sk 2 + s 2 ic 3 


c g (jq + sic 2 s 2 ic 3 ) 
s 3 + a 2 s 2 +a 2 s +a 0 


and the overall transfer function is 

yjs) _ g _ £o 

R{s) 1 + G H (s 3 + a 2 s 2 + a 2 s + a 0 ) + c 0 (Jc 3 s 3 + k 2 s + kj 

= 

s 3 + (a 2 + c 0 Jc 3 )s 2 + (a 3 + c 0 k 2 ) s + (a 0 + CqIcJ 


using the final value theorem. 


y BB (£) = lim 


lim y( t) : = y„ = lim s y(s) 

t -*• ff -0 

s c n J?(s) 


s-o s 3 + (a 2 + c 0 Jc 3 )s 2 + (a x + c 0 Jr 2 ) + (a 0 + c 0 k x ) 


for step input R(s) = s 1 
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•• y.s = 


a 0 + 


: = l for zero error 


ie, k^ — 1 ^ 


which is fixed. By appropriate selecting the value of k T , we 
implement any desired characteristic equation as we want. 


original plant: 


observer design 


dX 

dt 


= A X + h u 


y = £ t X 


Let X be the estimated state vector. 
Let 


can 
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=AX+hu+L ( y - 9 ) 

dt 

9 = q t 1 

where L is the observer matrix 

L - t h J 2 

define e. - X - X 

d_e dX _ d & 

dt dt dt 

= A X + h u - (A& + hu + L( £L t X - Q t &) ) 

= A( X - 1) - L £ r ( X - & ) 

= ( A - LQ T ) & 

By appropriate selecting the eigenvalues of e, the error of the 
estimated state vector will died out very quick. A state diagram is 
shown in figure 20. 


figure 20. 











Appendix IV 

Computer codes for nonlinear plant dyn amics 


{$N+,E+} 

program project (outputo) ; {3\21\1993} 

var 

d:char ; 

i, j, n, nstep, e, b: integer; 
h, h2 , Xlmax, X4max, t, r, u, dum:real? 
k : array [1. .3] of real; 
x, y : array [ 1 . . 2 , 1 . • 4 ] of real; 


f : array [ 1 . . 4 , 1 . . 4 ] of 
outputo : text ; 

{ This is a Runge-Kntta 
procedure initizing; 
begin 
t : =0 ; 

h:=l/1024 ; 
n : =4 ; 

for i:=l to 2 do 
for j:=l to 4 do 
begin 
x[i, j] :=0; 

y[i, j] :=o; 

end; 

for i:=l to 4 do 
for j:=l to 4 do 


real; 

method of order 4 . } 

{initial time} 
{incremental time step} 
{number of equations} 
{initial value} 
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f[i,j]:=o; 

nstep:=4096; {number of step} 

r:=l; 

k[l] :=1. 01; 
k[2] :=0. 021561; 
k[3] :=0. 0001505; 
end; 

procedure get_value_f; 
begin 

f [e, 1] := x [ b , 2 ] ; 

f [e, 2 ] := (382 . 399 6e- 6) * (4 . 3434e-3 + 

x[b,l]*x[b,3])*(0.00762*x[b,3] + 

0 . 57*x[b, 1] ) / (sqr(sqr( 0.007 62) - sqr (x[b, 1] ) ) ) ; 
dum := (y[b,2] - x[b,2])/h; {estimated the second derivate} 
u := r - (k[l]*x[b,l] + k[2]*x[b,2] + k[3]*dum); 

{control law} 

f [e, 3 ] := -44.2995*x[b,3] + 5.5374*u; 

f [e,4] := 1; { x[4]= time } 

end; 

procedure RK4SYS; 
begin 

h2:=0.5*h; 

for j:=0 to nstep do 
begin 
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for e:=l to 2 do 


begin 

b:=e; 

get_value_f ; 
for i:=l to n do 

x[2,i]:=x[l,i] + h2*f [e, i] ; 

end ; 
e : =3 ; 

get_value_f; {get f3> 
for i:=l to n do 

x[2,i]:=x[l,i] +h*f [e, i] ; 
e : =4 ; 

get_value_f; {get f4> 
y[l,2] :=x[l,2] ; 

Y [2 , 2 ] :=x[2 , 2 ] ; 

for i:=l to n-1 do {compute next x(t+h) } 

x[l,i]:=x[l,i]+h*(f[l,i]+2*(f[2,i]+f[3,i])+f[4,i])/6; 

x[l,4]:= t+j*h; {advance solution} 

if j mod 64 =0 then {write the result} 

begin 

for i:=l to n do 
write ( outputo ,x[l,i],','); 
writeln (outputo) ; 
writeln(x[l, 1] , ' , ' ,x[l,4] ) ; 
end{if loop} 
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end; {for j loop} 
end; {RK4SYS} 

begin {main} 

assign ( outputo , ' a : \pl . dat ' ) ; 

rewrite (outputo) ; 

writeln (outputo) ; 

initizing; 

repeat 

RK4SYS ; 

write('want change Y/N ?'); 
readln(d) ; 
if d= , y / then 
begin 

writeln ( ' kl= ' , k [ 1 ] , ' k2= / ,k[2],' k3=',k[3]); 
write (' enter kl,k2,k3 / ); 
readln (k[l] ,k[2] ,k[3]) ; 
end; 

until (do'y') ; 

writeln (outputo, ' Job completed. '); 
writeln (' Job completed.'); 
close (outputo) ; 
end. 
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